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(54) Abstract Trtle 

Signal analysis system and method for baseline optimisation 

(57) A method, apparatus, and article of manufacture for enhanced integration of a signal such as those 
generated by a chromatograph. In one embodiment an approximate baseline (1002) used for integration of a 
signal (101) is improved (1001) to reduce error. Alternate embodiments filter drift, noise, or both from a signal 
prior to integration or simulated distillation. The signal is processed according to a computer program stored 
in a memory. The methods include steps for subtracting (702) an approximate baseline (1002) from the signal 
(101), defining (703) a noise band (802-803) in the resulting difference (801), and forming a composite baseline 
(901) by substituting (700) signal data (101) for baseline data (1002) wherever the difference does not exceed 
the threshold (802-803), thus leaving the parts of the approximate baseline (1002) that correspond to siqnal 
peaks (102) unchanged (see eg. Fig 7). 
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SIGNAL ANALYSIS SYSTEM AND METHOD 

The present invention relates generally to signal analysis, especially two-dimensional 
signal analysis, and more particularly to a method, apparatus, and anicle of manufacture for 
the enhanced integration of time-varying signals which consist of one or more peaks rising 
above a background level containing noise and drift, such as those generated by 
chromatographic systems. 

A typical chromatographic analysis system includes an injection pon into which the 
sample is injected and mixed with an inert gas or liquid, a column through which the various 
dissolved components of the sample will travel at a rate related to the characteristics of the 
specific components, and a detector for measuring the concentration of each component as it 
exits the column and generating a time-varying signal corresponding to its concentration. The 
time between the injection of a sample and the detection of a specific component is called the 
retention time of that component. The time-varying signal is sent to a computing apparatus that 
integrates the signal by executing a computer program. The computing apparatus often provides 
a two-dimensional visual display of the signal, known in the art as a chromatogram The x-axis 
of the chromatogram represents time, and the y-axis typically corresponds to the amplitude of the 
signal 

A typical time-varying signal 101 generated by a gas chromatograph is shown in Figure 1. 
Peaks 102 in the signal correspond to the detection by the chromatograph of specific components 



and the amount of area under the peak and above the background level, correlate to the amount of 
the respective component present in the sample This information can be used for many 
applications, such as to check for the presence of, or to ensure the absence of. particular 
components in the sample, or to verify the proper concentration of the component in the sample. 

The computer program includes an integration method which detects peaks, determines 
the retention time 104 at which they occur, and determines their heights and areas To determine 
peak height and peak area, the computer must first determine a baseline for each peak. The 
baseline is a reference level that corresponds to the background level of the signaJ when no peaks 
are present The background level is typically a non-zero signal that may include noise (higher in 
frequency than the peaks) and drift (lower in frequency than the peaks; also known as wander). 

A chxomatographer visually drawing the baseline under a particular peak would 
extrapolate it through the center of the noise on both sides of the peak with a slope corresponding 
to the dnfi However, the presence of noise and drift in the background level make determining 
an accurate baseline difficult for a computer. As a result, prior an integration methods, in the 
presence of noise and drift, typically do not produce the same baseline as would a 
chxomatographer looking at the signal visually. 

A first pnor art integration method sets the baseline by looking at only a small region of 
the signal at a time As shown in Figure 2. this method sets the baseline 201 too low where noise 
202 is present because the method finds the lowest points near both s.des of the peak and draws 
the baseline 201 between these points. If. as shown in Figure 3. the detector produces a negative 
background signal disturbance 302 in the presence of two large peaks close together, this 
integration method will use the lowest point of the disturbance to set the baseline 301. Also. 
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misidentifying background noise as peaks, some valid peaks may not be detected This problem 
worsens as the signal-to-noise ratio decreases, as shown by the absence of any identified peaks 
for the signal 401 of Figure 4 

A second prior art method extends the baseline back in time from the average background 
level that occurs towards the end of an analysis in order to compensate for the effects of venting 
the solvent peak. As shown in Figure 5, this method sets the baseline 501 too low in the region 
where background drift 504 occurs after the venting of the solvent peak 502. In addition, if the 
background continues to drift throughout the run, this method may misinterpret a hump in the 
drift as a false peak 505. 

A third prior art method establishes a baseline using curve fitting techniques. However, 
the complexity of curve fitting increases with the complexity of the baseline shape, which makes 
such a method difficult to implement and requires excessive computational power when the 
signal contains noise and drift. To reduce complexity, this method may alternatively divide the 
baseline into multiple segments to allow simpler, piecewise curve fitting. However, this 
modification introduces discontinuities between curve segments that affect the accuracy of 
integration results. 

A fourth prior an method filters noise by performing a fourier transform However, it is 
difficult to filter out all the noise without distorting the peaks, since typically there is overlap in 
frequency content between the peaks and the noise 

A chromatogram which displays an accurate baseline along with the signal is readily 
usable by an expert chromatographer. However, chromatograms are sometimes analyzed by non- 
chromatographers, for whom baseline drift and noise can cause confusion A better display for 



filtered, leaving only undistorted peaks and a flattened baseline. A filtered chromatogram could 

be accurately integrated using simple integration methods, and would provide more accurate 

results when analyzed using simulated distillation techniques. 

In the analysis of signals consisting of one or more peaks rising above a background level 

which contains noise, drift, or a combination of both, there has been a need for a simple method 
of determining an accurate baseline which eliminates noise and drift without distorting the peaks, 
in order to calculate accurate peak heights and areas. There has aJso been a corresponding need 
for filtering the signal to eliminate noise, drift, or both, prior to display, integration, or simulated 
distillation. 



The present invention may be implemented as a method of optimizing the baseline of a 
signal to account for noise and drift without introducing peak distortion. This method facilitates 
accurate measurements of peak height and peak area during integration, without undue 
complexity The present invention may alternatively be implemented as a simple method of 
filtering the signal to attenuate noise, drift, or both A chomatogram filtered according to this 
invention communicates information more clearly to non-chromatographers. can be accurately 
integrated by simple integration methods, and prov.des accurate results in simulated distillation 
analysis The invention may be embodied in a program storage medium readable by a computer 
which contains program code that performs the baseline optimization and signal filtering 
methods It may also be embodied in a computing apparatus having a computer program which 
performs the baseline optimization and signal filtering methods; a gas chromatographic apparatus 




A preferred first embodiment of the present invention is a method that produces an 
improved baseline from an approximate baseline by subtracting the approximate baseline from 
the signal, defining a noise band in the resulting difference, forming a composite baseline by 
substituting signal data for approximate baseline data at those segments of the baseline where the 
difference falls within the noise band (thus leaving the segments of the approximate baseline that 
correspond to the signal peaks unchanged), and smoothing the composite baseline to form the 
improved baseline. These steps may optionally be performed repetitively to further improve the 
baseline, using the improved baseline from the prior iteration as the approximate baseline for the 
subsequent iteration. A further option is to integrate the signal using the improved baseline to 
obtain the peak height and peak area. 

A preferred second embodiment of the present invention is a method that filters a signal 
to attenuate drift by smoothing the signal to form an approximate baseline, subtracting the 
approximate baseline from the signal, defining a noise band in the resulting difference, forming a 
composite baseline by substituting signal data for approximate baseline data at those segments of 
the baseline where the difference falls within the noise band (thus leaving the segments of the 
approximate baseline that correspond to the signal peaks unchanged), smoothing the composite 
baseline, and subtracting the smoothed composite baseline from the initial signal to form the 
filtered signal with drift attenuated These steps may optionally be performed repetitively to 
further filter the drift, using the smoothed composite baseline from the prior iteration as the 
approximate baseline for the subsequent aeration. A further option is to integrate the filtered 
signal to obtain the peak height and peak area once the baseline has been improved using a 
simple integration algorithm. 
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A preferred third embodiment of the present invention is a method that filters a signal to 
^ attenuate both noise and drift by smoothing the signal to form an approximate baseline, 
subtracting the approximate baseline from the signal, defining a noise band in the resulting 
difference, forming a composite baseline by substituting signal data for approximate baseline 
data at those segments of the baseline where the difference falls within the noise band (thus 
leaving the parts of the baseline that correspond to the signal peaks unchanged), and subtracting 
the composite baseline from the initial signal to form the filtered signal with noise and drift 
attenuated. These steps may optionally be performed repetitively to further filter the noise and 
drift, by smoothing the composite baseline from the prior iteration and then using it as the 
approximate baseline for the subsequent iteration. A further option is to integrate the filtered 
signal to obtain the peak height and peak area once the baseline has been improved using a 
simple integration algorithm. Alternatively, a simulated distillation analysis may be performed 
on the filtered signal. 

A preferred fourth embodiment of the present invention is a method that filters a signal to 
attenuate noise by smoothing the signal to form an approximate baseline, subtracting the 
approximate baseline from the signal, defining a noise band in the resulting difference, forming 
an unsmoothed composite baseline by substituting signal data for approximate baseline data at 
those segments of the baseline where the difference falls within the no.se band (thus leaves the 
parts of the baseline that correspond to the signal peaks unchanged), smoothing the composite 
baseline to form a smoothed composite baseline, and subtracting the unsmoothed composite 
baseline from the initial signal while adding the smoothed composite baseline to the initial signal 
to form the filtered signal with noise attenuated. These steps may optionally be performed 




iteration as the approximate baseline for the subsequent iteration. 

Other aspects and advantages of the present invention will become apparent from the 
following detailed description, taken in conjunction with the accompanying drawings, illustrating 
by way of example the principles of the invention. The claims alone, not the preceding summary 
or the following detailed description, define the invention. 

Fig. 1 is an exemplary chromatogram showing a chromatographic signal including peaks 
and noise, and an approximate baseline. 

Fig. 2 is an exemplary chromatogram showing a chromatographic signal having 
background noise and an inaccurate baseline resulting from use of a prior art integration method. 

Fig. 3 is an exemplary chromatogram showing a chromatographic signal having a 
negative background disturbance and an inaccurate baseline resulting from use of a prior an 
integration method. 

Fig. 4 is an exemplary chromatogram showing a chromatographic signal with a low 

signal-to-noise ratio where no peaks are detected by a prior an integration method. 

Fig. 5 is an exemplary chromatogram showing a chromatographic signal including peaks, 

noise, and drift. The chromatogram shows an inaccurate baseline and a false peak resulting from 

use of a prior an integration method- 
Fig. 6 is a block diagram of a gas chromatographic instrumentation system suitable for 

the practice of the present invention. 

Fig. 7 is a simplified representation of an algorithm for performing baseline optimization 




Fig. 8 is an exemplary chromatogram showing the baseline-adjusted signal generated b> 
subtracting the approximate baseline from the signal of Fig I, and showing the threshold values 
defining the noise band for the signal. 

Fig 9 is an exemplary- chromatogram showing the composite baseline (before smoothing) 
generated by the first iteration of the baseline optimization algorithm operating on the signal of 
Fig 1. 

Fig. 10 is an exemplary chromatogram showing, on an expanded scale relative to Fig 1, 
the initial approximate baseline, the first pass improved baseline, and the second pass improved 
baseline generated by the baseline optimization algorithm operating on the signal of Fig. 1. 

Fig 1 1 is an exemplary chromatogram showing the baseline-adjusted signal and the 
threshold values defining the noise band for the signal of Fig. 1, after the second iteration of the 
baseline optimization algorithm. 

Fig 12 is a simplified representation of an algorithm for performing signal filtering to 
attenuate drift in accordance with the present invention. 

Fig 13 is an exemplary chromatogram showing the filtered signal resulting from applying 
the signal filtering method of Fig. 12 to the signal of Fig. 5. 

Fig 14 is a simplified representation of an algorithm for performing signal filtering to 
attenuate drift and noise in accordance with the present invention. 

Fig 15 is an exemplary chromatogram showing the composite baseline (offset slightly 
downward for viewability ) generated by applying the signal filtering method of Fig. 14 to the 
signal of Fig 5. 

Fig. 16 is an exemplary chromatogram showing the filtered signal resulting from applying 



Fig 17 is a simplified representation of an algorithm for performing signal filtering to 
attenuate noise in accordance with the present invention. 

Fig. 18 is an exemplary chromatogram showing the filtered signal resulting from applyin 
the signal filtering method of Fig. 17 to the signal of Fig. 5 

As shown in the drawings for purposes of illustration and as described in the following 
detailed description, the method, apparatus, and article of manufacture of the present invention 
improve the analysis of electronic signals consisting of one or more peaks greater in amplitude 
than a background level containing noise higher in frequency than the peaks, drift lower in 
frequency that the peaks, or both noise and drift. According to the present invention, integration 
and simulated distillation of these signals produce more accurate measurements, and the 
chromatograms of the analyzed signals can be understood more readily by non-experts. Prior 
methods of analyzing these signals have been unduly complex or have yielded less accurate 
measurements 

Electronic signals usable with the present invention can be generated by a variety of 
instrumentation, and may represent a variety of physical, electrical, or chemical properties. 
Examples of such instrumentation include, but are not limited to, gas and liquid chromatographs, 
spectrum analyzers, oscilloscopes, and seismographs Signals from some instruments may be 
generally bipolar, having sinusoidally-shaped peaks centered around the background level. The 
signals from chromatographs are generally unipolar, having Gaussian-shaped peaks rising above 
the background level (however, signals dropping below the background level, often called 
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signals according to the preferred practice of the invention, and therefore the deta.led description 
of the invention will be directed to a gas chromatographic analytical system However, it should 
be understood that the teachings herein are applicable to any electronic signals having the 
amplitude and frequency properties described above. 

Gas chromatographs usable in the practice of the present invention include the Hewlett- 
Packard 6890 series and 5890 series, while liquid chromatographs include the Hewlett-Packard 
1 100 senes and 1050 series. 

System, Apparatus, and Program Storage Medium 

The preferred embodiment of an analytical instrument system 600 illustrated in Figure 6 
includes a gas chromatograph 610, and a computer system 620 having a signaJ input 621, 
processor 622. memory 623, interface 624. input device 627. and output device 628. In order to 
perform a chromatographic separation of a sample 601. a quantity of the sample 601 is injected 
into a fluid stream, preferably in the form of a pressurized carrier gas. via an injection port 61 1 
supplied with carrier gas stream 613. The injection port 61 1 provides a portion of the 
sample/carrier gas mixture 612 to a separation column 614. The column 614 ,s positioned within 
a temperature-controlled thermal chamber or oven 615. The carrier gas/sample combination 612 
passing through the column 614 is exposed to a temperature profile resulting from the operation 
of the heater 616 w ithin the oven 614 

As the carrier gas (containing the separated components of the sample) exits the column 
614. the presence of one or more of the constituent components of the sample is detected by a 
detector 6 1 7. The detector 61 7 can be any of the GC detectors known in the an. as long as it is 
capable of determining at least one physicocherrucal property of the earner gas which exits the 




Such GC detectors include the flame ionization detector (FID), photoionization detector 
(PID), nitrogen phosphorous detector (NPD), flame photometric detector (FPD), thermal 
conductivity detector (TCD), atomic emission detector (AED), electrolytic conductivity detector 
(ELCD), and electron capture detector (ECD). Mass spectral detectors and infrared spectral 
detectors are also known. 

The detector output signal 618 may be either in analog or digital electronic form The 
detector output signal is then received by the signal input 621 for communication to the processor 
622 If ihe signal is in analog form, the signal input may include an analog-to-digital converter 
(not shown) to convert the signal to a computer-readable form. The converted detector output 
signal is typically stored in the memory 623. Preferably, the detector output signal is provided in 
the form of data representative of one or more peaks each corresponding to detection of a 
constituent component of the sample, the peak(s) rising above a background signal level 
containing noise, drift, or both noise and drift. The operations performed by the computer system 
620 on the signal data in accordance with the present invention will be described subsequently 

The processor 622 may include computing devices amenable to the practice of this 
invention, e.g , one or more computing devices such as computers, microprocessors, 
microcontrollers, switches, logic gates, or any equivalent logic device capable of performing the 
computations described herein below. The processor 622 preferably is coupled with the interface 
624, input device 627, and output device 628 Input devices preferably include one or more of a 
keyboard, keypad, pointing device, or remote computer (not shown) for inputting operating 
condition parameters, system calibration data, and the like. Output devices preferably include 
alphanumeric or video displays, printers, or a remote computer (not shown). 
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The preferred computer system 620 further includes memory 623 in the form of volatile 
and non-volatile devices in which computer programs 626, input and output information, 
operating condition parameters, and system information can be stored and retrieved 
Contemplated computer programs 626 include method embodiments of the present invention 
Contemplated memory devices include, but are not limited to, ROM, RAM, floppy disk, hard 
disk, and CD-ROM. Operating commands, detector response attributes, column temperature 
programs, and other information necessary to perform a chromatographic analysis may be entered 
into the processor 622 by way of an input device 627 or retrieved from memory 623. 

The chromatogram, peak heights, and peak areas calculated by the processor 622 for the 
signal may be transmitted to an output device 628 and displayed or printed. The interface 624 
may further comprise network and bus system (input/output or I/O) controllers, isolation devices, 
clocks, and other related electronic components for performing control, processing, and 
communication tasks other than those described herein. A remote computing apparatus (not 
shown) may connect to the interface 624 and cooperate with the computer system 620 to 
implement the present invention as a client/server application. 

The method embodiments of the present invention are operable in the system 600, and 
may be contained in the memory 623. 

Signal Data Representation and Storage 

The signal data, herein referred to mathematical^ as s(x). is stored in the memory 623 in 
a format usable for two-dimensional data. The dependent variable s(x) represents the signal in a 
chromatogram, and the independent variable x represents time. In a first embodiment where the 
values of x represent a fixed interval of time, s(x) may be stored in a data structure comprising a 




represent a variable interval of time, a data structure comprising an array of ordered pairs (x.s(x)) 
may be used. Data structures of the corresponding type may be used to contain the other 
dependent variable data disclosed in the remainder of the detailed description, such as baselines 
b(x) and baseline-adjusted signals a(x). 

Smoothing Algorithm 

Several embodiments of the present invention include a method step for smoothing a 
signal or a baseline. The purpose of this step is to attenuate high frequencies (noise and peaks) in 
the signaJ or baseline while having little or no effect on low frequencies (drift). Smoothing 
conserves area; when a peak gets smoothed, its area gets redistributed over a wider range. 
Visually, the peak becomes lower in amplitude but extends over a wider period of time. 

A variety of high-frequency filtering algorithms may be utilized to perform the smoothing 
step The general approach to smoothing is to replace each of the original data points with the 
weighted average of it and its nearest neighbors. The number of points averaged, and the 
weighting factors applied to each point t determine the frequency-response curve. There is one 
complication to any smoothing algorithm: data points near the start and end of the data do not 
have the required number of nearest neighbors on one side, so their smoothed values cannot be 
calculated in the usual way. Often the solution to this problem is to have fewer points in the 
smoothed data than were in the original data. The first and last (n-l)/2 points (where n is the 
number of points averaged) are simply missing Where multi-pass smoothing is applied, more 
data will being lost at each subsequent pass. 

The preferred smoothing algorithm uses a simple "boxcar" filter in which all of the 
weighting factors are the same and the only variable is the number of points averaged. To 




points averaged. This achieves high-frequency rejection with less distortion, and is funciionallv 
equivalent to smoothing with unequal weighting factors To avoid data loss at the start and end 
of the signal or baseline, the algorithm extrapolates additional data points prior to smoothing 
which axe included in the filtering. 

An alternate smoothing embodiment having unequal weighting factors uses aSavitzkv- 
Golay (least-squares) filter. The variables in least-squares smoothing are the number of points 
averaged and the order of equation to which the data is to be fit. Although least-squares 
smoothing is more complex than boxcar smoothing, it can achieve a much better tradeoff 
between high-frequency attenuation and data distortion in a single pass. Data loss at the start and 
end of the signal or baseline can be avoided by using weighting factors that calculate the least- 
squares value at a point other than the center point for each data point near the ends of the data. 

Where an embodiment of the present invention calls for smoothing a signal that contains 
large, broad peaks to form an initial or approximate baseline, a more accurate initial baseline may 
be constructed after smoothing by forming a composite consisting of smoothed signal at those 
values of x where the signal approximately represents the background, and straight line segments 
between the endpoints of the smoothed signal to bridge the gaps where the signal approximately 
represents peaks. In order to determine where the signal represents peaks, the steps of forming a 
baseline-adjusted signal 702 and measuring the noise to set threshold values 703 may be 
performed, as shown in Figure 7 and described subsequently. 

Noise Measuring Algorithm 
Several embodiments of the present invention include a method step for measuring the 
noise in a baseline-adjusted signal to define positive and negative threshold values, t P and t N . 
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contiguous groups of data points and measuring the noise in each segment In each segment, 
noise can be measured as either peak-to-peak or RMS noise Peak-to-peak noise is the difference 
between the largest and the smallest values within the segment RMS is the standard deviation of 
all of the data values within the group. Peak-to-peak is faster and easier to measure, but RMS is 
more accurate Peak-to-peak noise is roughly six times larger than RMS noise for the same data. 
The preferred embodiment measures RMS noise. The thresholds t P and t s are set to plus and 
minus a factor times the RMS noise respectively. The preferred embodiment uses a factor of 2.0. 

In measuring noise, the algorithm must ensure that it does not confuse peaks or drift with 
noise In the preferred embodiment of noise measurement, any drift in the input signal has 
previously been subtracted from the baseline-adjusted signal in which the noise is measured. In 
order to exclude peaks, the preferred embodiment finds those segments with the smallest non- 
zero noise (background noise is assumed to be much lower than peaks). Once the noise 
measurements have been made for all segments, the smallest value can be chosen as the baseline 
noise Since a particular detector may saturate above a certain signal level, causing it to have 
erroneously low noise, the preferred embodiment rejects any measured noise value of zero and 
selects the smallest non-zero value as the noise. 

A further aspect of the noise measurement algorithm involves which groups of points are 
used to perform the measurement. The grouping can be as simple as dividing up the data points 
into some number of approximately equal groups (for example, one hundred groups of n/100 
each, where n is the number of points in the signal) and calculating the noise of each group. A 
more sophisticated approach measures noise on more than 100 groups but still includes 1/100 of 
the data points in each group, for example the first group might be points 1 to 100. the second 
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increments of the starting data point (for example, points l to 1 00, 2 to 1 0 1 , 3 to 1 02, etc.). 

Baseline Optimization Method 
A preferred first embodiment of the present invention is a novel method of optimizing a 
baseline used for integration of a signal. The general approach is to begin with an initial 
approximate baseline that approximately tracks the background drift of the signaJ but has an error 
component, and improve the accuracy of this baseline by forming a composite baseline which 
tracks the signaJ noise and then smoothing the composite baseline to center the baseline within 
the noise. The error component of this improved baseline is smaller than that of the initial 
baseline. 

As illustrated by applying the method of Figure 7 to the signal data of Figure l by way of 
example, the first embodiment begins with a initial step 701 in which an approximate baseline 
b(x) 105 is generated- The approximate baseline may be generated by smoothing the signal s(x), 
or by other methods which produce a baseline that generally tracks the background drift of the 
signal. The approximate baseline contains inaccuracies, particularly in the region of signal peaks 
and negative background disturbances, due to the area conservation of the smoothing algorithm 

Once the approximate baseline has been obtained, the following optimization steps can be 
performed one or more times to optimize the baseline. The first step 702 in optimization is to 
subtract the baseline from the signal, the difference forming a baseline-adjusted signaJ a(x) 801 
centered around zero counts, as shown in Figure 8. The baseline-adjusted signal shows the 
deviation of the signal from the baseline, for example, at a value of x where the baseline b(x) 
equals the signal s(x), the base line -adjusted signal a(x) equals zero counts. Since the smoothed 
baseline from the prior step tracked the drift, drift will be attenuated in the baseline-adjusted 




The second optimizing step 703 js to calculate the noise in the baseline-adjusted signal 
a(x) 801 and set positive (t P ) 802 and negative (t s ) 803 threshold values outside the range of, and 
surrounding, the noise. Signal data s(x) at those values of x where the baseline-adjusted signal 
falls within the noise band defined by the thresholds approximately represent signal background. 
Signal data at those values of x where the baseline-adjusted signal fails outside this band 
approximately represent signal peaks or negative background disturbances. 

The third step 704 is to substitute data from the signal s(x) into the baseline b(xj at those 
values of x where the baseline-adjusted signal a(x) falls within the noise band defined by the 
threshold values. After this step has been performed, the baseline 901 in Figure 9 becomes a 
composite of signal data 902 at those values of x where the signal approximately consists of 
background, and previous baseline 903 at those values of x where the signal approximately 
consists of peaks or negative background disturbances. As will be appreciated by comparing 
Figure 9 with Figure 1, the composite baseline 901 approximates signal 101 data, with peaks 102 
and disturbances 106 outside the noise band clipped. Replacing baseline data b(x) with signal 
data s(x) has the effect of removing from the composite baseline 901 much of the area of the 
peaks which was included in the baseline 105 following smoothing 701. 

The fourth step 705 is to smooth the composite baseline 901 to form the improved 
baseline 1001 which has a reduced error component. Since signal data representing peaks and 
disturbances outside the threshold values have been removed from the composite baseline 901, 
they do not distort the first pass improved baseline 1001 in Figure 10 resulting from the 
smoothing step 705. It will be appreciated that the first pass improved baseline 1001 is closer to 
a baseline centered through the noise and drift than the initial baseline 1002, and is closer to a 
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The preceding opt.mizing steps 702. 703. 705. and 705 may be performed iteratively 706 
to further improve the baseline, substituting the improved baseline generated by the previous 
iteration as the approximate baseline for the next iteration. Because the first pass improved 
baseline 1001 is more accurate than the initial baseline 1002. the deviations from zero counts in 
the second pass baseline-adjusted signal 1 101 of Figure 1 1 are further reduced, allowing second- 
pass thresholds t P 1 102 and t N 1103 to be set tighter relative to first-pass thresholds t P 802 and t N 
803. which in turn results in a more accurate second-pass composite baseline. The second 
iteration yields the second pass improved baseline 1003 of Figure 10. For a fairly simple signal 
with small peaks and disturbances, good results are possible after only one or two iterations, but 
to handle the more generaJ case the preferred embodiment iterates eight times. Additional 
iterations of the optimization steps yield further improvement, at the cost of increased processing 



time. 



use 



After the baseline has been optimized, the signal may be integrated via a prior an 
integration algorithm that uses this optimized baseline to define the time and signal level to 
for the starting and ending baseline points for each peak so as to more accurately calculate peak 
heights and areas 

Signal Filtering Method to Attenuate Drift 

As illustrated by applying the method of Figure 12 to the signal data of Figure 5 by way 
of example, a preferred second embodiment of the present invention is a novel method of 
filtering a signal 503 to attenuate background drift 504. The general approach used in the second 
embodiment is to generate a smoothed composite baseline which tracks the signal drift, and then 
subtract the baseline from the signal to attenuate the drift The sirru laxities between this method 
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This embodiment begins with a step 1201 in which an initial baseline bfx) is generated by 
smoothing the signal s(x) 503. The smoothing operation generates a baseline that approximately 
tracks low-frequency phenomena such as background drift but not higher-frequency phenomena 
such as background noise and peaks. 

The second step 1202 is to subtract the baseline b(x) from the signal s(x) 503, the 
difference forming a baseline-adjusted signal a(x) centered around zero counts. Since the 
smoothed baseline b(x) tracks drift in the signal but not peaks or noise, it will be appreciated that 
drift in the baseline-adjusted signal a(x) is attenuated while peaks and noise remain. 

The third step 1203 is to calculate the noise in the baseline-adjusted signal a(x) and set 
positive (tp) and negative (t N ) threshold values outside the range of, and surrounding, the noise. 
Signal data at those values of x where the baseline-adjusted signal falls within the noise band 
defined by the thresholds approximately represent background. Signal data at those values of x 
where the baseline-adjusted signal falls outside this band approximately represent peaks. 

The fourth step 1204 is to substitute data from the signal s(x) into the baseline b(x) at 
values of x where the baseline-adjusted signal a(x) falls within the noise band defined by the 
threshold values. Following the performance of this step, the baseline becomes a composite of 
signal data at those values of x where the signal approximately consists of background, and 
previous baseline at those values of x where the signal approximately consists of peaks The 
composite baseline approximates signal data with any peaks outside the noise band clipped. It 
will be appreciated that substituting signal data into the baseline makes the composite baseline 
approximately track the drift and noise in the signal. 
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The fifth step 1205 is to smooth the composite baseline b(x). The smoothing operation 
results in a baseline b(x) that still approximately tracks background drift, but the effects of noise 
get smoothed out and attenuated. 

The second, third, fourth, and fifth steps 1202. 1203. 1204, and 1205 above may be 
performed iteratively 1207 to improve how well the baseline b(x) tracks the background drift, by 
using the smoothed composite baseline resulting from the previous iteration as the approximate 
baseline for the next iteration. For a fairly simple signal with small peaks and disturbances, good 
results axe possible after only one or two iterations, but to handle the more general case the 
preferred embodiment iterates eight times. Additional iterations further attenuate drift, at the 
cost of increased processing time. 

After the final iteration has been performed, the last step 1206 is to subtract the baseline 
b(x) from the signal s(x) to form the filtered signal 1301 shown in Figure 13. It will be 
appreciated that this results in a signal 1301. centered around zero counts, in which background 
drift has been attenuated but background noise has not. The heights of the peaks 1302 have been 
adjusted to compensate for drift, as can be seen by comparison to the peaks 506. Peak area is 
conserved 

The filtered signal may further be displayed as a chromatogram, or u may be integrated 
via a prior an integration algorithm to calculate peak heights and areas. Since drift has been 
filtered from the signal, integration is simplified in that a signaJ level of zero counts may be used 
for the starting and ending baseline points for ail peaks In addition, false peaks due to drift are 
eliminated, as can be seen by theabsense in Figure 13 of the false peak 505. 
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Signal Filtering Method to Attenuate Noise and Drift 

As illustrated by applying the method of Figure 14 to the signal data of Figure 5 by way 
of example, a preferred third embodiment of the present invention is a novel method of filtering a 
signal 503 to attenuate both noise 507 and drift 504 The general approach is to generate a 
smoothed baseline which tracks the signal drift, then create a composite baseline which tracks the 
signal noise, and finally subtract the composite baseline from the signal to attenuate the drift and 
noise. The similarities between this method and the signal filtering method to attenuate drift 
described previously will be appreciated by those skilled in the an. 

This embodiment begins with a step 1401 in which an initial approximate baseline b(x) is 
generated by smoothing the signal s(x) 503. The smoothing operation generates a baseline that 
tracks low-frequency phenomena such as drift but not higher-frequency phenomena such as noise 
and peaks. 

The second step 1402 is to subtract the baseline b(x) from the signal s(x) 503, the 
difference forming a baseline-adjusted signal a(x) centered around zero counts. Since the 
smoothed baseline b(x) tracks drift in the signal but not peaks or noise, it will be appreciated that 
drift in the baseline-adjusted signal a(x) is attenuated while peaks and noise remain 

The third step 1403 is to calculate the noise in the baseline-adjusted signal a(x) and set 
positive (t P ) and negative (t N ) threshold values outside the range of. and surrounding, the noise 
Signal data at those values of x where the baseline-adjusted signal falls w ithin the noise band 
defined by the thresholds approximately represent background. Signal data at those values of x 
where the baseline-adjusted signal falls outside this band approximately represent peaks. 

The fourth step 1404 is to substitute data from the signal s(x) into the baseline b(x) at 
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threshold values. Follow.ng the performance of th.s step, the baseline becomes a composite of 
s.gnal data at those values of x where the s.gnal approx.mately cons.sts of background, and 
previous baseline at those values of x where the s.gnal approx.mately cons.sts of peaks. The 
composite basel.ne approximates s.gnal data with any peaks outs.de the no.se band clipped. It 
will be appreciated that substituting s.gnal data into the baseline makes the composite basel.ne 
approximately track the drift and noise in the signal. 

The second, third, and fourth steps above may be performed iterat.vely 1406 to further 
improve how well the compostte basel.ne tracks the background drift and noise. If the method ,s 
performed iterat.vely. the compostte baseline b(x) generated by the fourth step 1404 of rh, nrinr 

- I *" 

iteration will be smoothed 1405 and then will replace the approximate basel.ne used for 
performing the subtraction step 1402 of the next iteration For a fairly simple signal with small 
peaks and disturbances, good results are possible after only one or two aerations, but to handle 
the more generaJ case the preferred embodiment iterates eight times AdditionaJ iterations 
further attenuate drift and no.se. at the cost of increased processing time. Figure 15 shows the 
composite basel.ne 1501 produced by the final iteration/overlaid with the s.gnal 1502 (for 
purposes of clarity, the basel.ne is illustrated offset slightly downward) It w,H be appreciated 
that the basel.ne closely tracks the s.gnal everywhere but where peaks are present. 

Once iterat.on has been completed and the final composite baseline b(x) has been 
obtained, the compos.te basel.ne ,s subtracted 1407 from the s.gnal s(x). As shown in F.gure 16, 
the background drift and no.se are attenuated, m fact, after sufficient aerations the resulung 
signal s(x) 1601 will be zero counts everywhere except where peaks 1602 exist. The height of 
the peaks 1602 above the basel.ne is unchanged, as can be seen by comparison to the peaks 506. 
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and peak area is conserved. In addition, false peaks due to drift are eliminated, as can be seen by 
the absence in Figure 16 of the false peak 505 

The drift-adjusted signal may further be displayed as a simple chromatogram, or it may be 
integrated via a prior art integration algorithm to calculate peak heights and areas. Since drift 
and noise have been filtered from the signal, integration is simplified. A signal level of zero 
counts is used for the starting and ending baseline points of all peaks, and a simple peak detector 
will suffice since any positive deviation from zero counts indicates a peak. In addition, false 
peaks due to drift are eliminated, as can be seen by the absence in Figure 16 of the false peak 
505 

Another application of this filtering method is in simulated distillation, which does not 
involve integration. In simulated distillation the signal is divided up into small segments called 
area slices Each area slice corresponds to a certain boiling range of a hydrocarbon mixture. 
Which peaks make up the area within a slice is not important, since all peaks in the slice have 
approximately the same boiling point; rather, it is the total area under the peaks of each slice that 
are measured. Significant errors result from simulated distillation when noise and drift are 
significant relative to the signal. However, by using this filtering method to remove drift and 
noise from the signal prior to performing simulated distillation, the area slices will only contain 
the area of the peaks This technique is only applicable where there are discrete peaks, as in 
lower-boiling mixtures Higher boiling ranges tend to produce "humps" of highly merged peaks 
and the filtering method will treat the humps as drift and subtract out much of their area 

Signal Filtering to Attenuate Noise 
As illustrated by applying the method of Figure 17 to the signal data of Figure 5 by way 
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a signal 503 to attenuate noise 507. The general approach ts to form a composite baseline which 
tracks noise and dnft. then smooth the composite baseline to form an intermediate baseline 
which tracks only drift, subtract the composite baseline from the signaJ to attenuate both noise 
and drift, and finally add the intermediate baseline back in to reconstitute the drift Those skilled 
in the art will appreciate that this method combines elements of the two signal filtering methods 
described above. 

This embodiment begins with a step 1701 in which an initial baseline b(x) is generated by 
smoothing the signaJ s(x) 503. The smoothing operation generates a baseline that tracks low- 
frequency phenomena such as drift but not higher-frequency phenomena such as noise and peaks. 

The second step 1702 is to subtract the baseline b(x) from the signal s(x). the difference 
forming a baseline-adjusted signaJ a(x) centered around zero counts Since the smoothed 
baseline b(x) tracks drift in the signal but not peaks or noise, it will be appreciated that the drift 
in the baseline-adjusted signaJ a(x) is attenuated while peaks and noise remain. 

The third step 1703 is to calculate the noise in the baseline-adjusted signal a(x) and set 
positive (t P ) and negative (t N ) threshold values outside the range of, and surrounding, the noise. 
Signal data at those values of x where the baseline-adjusted signal fails within the noise band 
defined by the thresholds approximately represent background. Signal data at those values of x 
where the baseline-adjusted signal falls outside this band approximately represent peaks. 

The fourth step 1704 is to substitute data from the signal s(x) into the baseline b(x) at 
values of x where the baseline-adjusted signal a(x) falls within the noise band defined by the 
threshold values. Following the performance of this step, the baseline becomes a composite of 
signal data at those values of x where the signal approximately consists of background, and 




composite baseline approximates signal data with any peaks outside the noise band clipped It 
will be appreciated that substituting signal data into the baseline makes the composite baseline 
approximately track the background drift and noise in the signal 

The fifth step 1705 is to smooth the composite baseline b(x) to form an intermediate 
baseline b2(x) that still approximately tracks background drift, but no longer tracks background 
noise due to effects of the smoothing operation. 

The second, third, fourth, and fifth steps above may be performed iteratively 1706 to 
further improve how well the composite baseline b(x) and the intermediate baseline b2(x) track 
the background. If the method is performed iteratively, the intermediate baseline b2(x) generated 
by the fifth step 1705 of the prior iteration will replace the baseline b(x) 1707 for the subtraction 
step 1702 of the next iteration. For a fairly simple signal with small peaks and disturbances, 
good results are possible after only one or two iterations, but to handle the more general case the 
preferred embodiment iterates eight times. Additional iterations further attenuate noise, at the 
cost of increased processing time. 

Once iteration has been completed and the final composite baseline b(x) and intermediate 
baseline b2(x) have been obtained, the composite baseline is subtracted 1708 from the signal s(x) 
to attenuate the background drift and noise, and the intermediate baseline b2(x) is added 1708 to 
the signal s(x) to add the background drift back in As shown in Figure 18, the background noise 
is attenuated while the drift 1 803 is generally unaffected As a result, the amplitude of peaks 
1802 is substantially the same as the peaks 506 in the original signal 503. It will be appreciated 
that because of the addition of the background drift, the resulting signal s(x) 1801 will not be 
centered around zero counts, but around the level of the original background drift 1803. 



26 

The drift-adjusted signal may further be displayed as a chromatogram which show;, the 
drift in the signal without the masking effects of noise The signal may also be integrated via a 
prior art integration algorithm to calculate peak heights and areas 

Conclusion 

From the foregoing it will be apparent that the invention provides novel and 
advantageous methods of optimizing the baseline of a signal to properly account for the 
effects of noise and drift, and of Filtering drift, noise, or both from a signal prior to 
performing subsequent processing operations. An accurate baseline ensures proper peak 
recognition and produces accurate integration results for peak height and peak area Pre- 
fihenng the signal to remove drift and noise provides a visually cleaner display, allows 
integration to be performed by simpler and faster algorithms, and allows accurate 
simulated distillation to be performed. 

Although several specific embodiments of the invention have been described and 
illustrated in the area of gas chromatography, the invention is not to be limited to the 
specific forms, methods, or arrangements of pans so described and illustrated For 
example, the embodiments of the invention described herein may operate upon signals 
which are generated by other electronic instruments such as spectrum analyzers and 
seismographs. The invention is defined by the claims. 




CLAIMS 

1 . A method for analyzing an input signal having a background level, comprising the method 
steps of: 

acquiring the input signal (101); 

subtracting (702) an approximate baseline (1002) having a first error component from the 
input signal (101) to form a baseline-adjusted signaJ (801) having baseline -adjusted noise; 

calculating (703) first (802) and second (803) threshold values surrounding the range of 
the baseline-adjusted noise and defining a noise band between the first and second threshold 
values; and 
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baseline-adjusted signal (801 ) falls within the noise band to form a composite baseline (901). 

2 The method for analyzing an input signal as claimed in claim 1, further comprising the method 
step of: 

smoothing (701) the input signal (101) to form the approximate baseline (1002). 

3. The method for analyzing an input signal (503) as claimed in claim 1, wherein the background 
level includes a first drift component, further comprising the method steps of: 

smoothing (1205) the composite baseline to form an improved baseline; and 
subtracting (1206) the improved baseline from the input signal to form a filtered signal 
(1301 ) having a second drift component attenuated relative to the first drift component. 

4. The method for analyzing an input signal as claimed in claim 3, further comprising the 
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method steps of subtracting the approximate baseline, caJcuJating first and second threshold 
values, substituting the input s.gnai into the approximate baseline, smoothing the composite 
baseline, and replacing the approximate baseline are performed iteratively (1207), whereby the 
second drift component is further attenuated with each iteration. 

5. The method for analyzing an input signal (503) as claimed in claim 1. wherein the background 
level includes a first drift component and a first noise component, further comprising the method 
step of subtracting (1407) the composite baseline from the input signal to form a filtered signal 
(1 60 1 ) having a second drift component attenuated relative m rfv fi~t . _ , 
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second noise component attenuated relative to the first noise component. 

6. The method for analyzing an input signal as claimed in claim 5, further comprising the 
method steps of: 

smoothing (1405) the composite baseline to form an improved baseline, and 
replacing the approximate baseline with the improved baseline; 

wherein the method steps of subtracting the approximate baseline, calculating first and 
second threshold values, substituting the input signal into the approximate baseline, smoothing 
the composite baseline, and replacing the approximate baseline are performed iteratively (1406), 
whereby the second drift component and second no.se component are further attenuated with 

each iteration. 
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7 The method for analyzing an input signal as claimed in claim 5, wherein the input signal is 
chromatographic data, further comprising the method step of performing simulated distillation 
calculations on the filtered signal. 



8. The method for analyzing an input signal (503) as claimed in claim 1, wherein the background 

level includes a first noise component, further comprising the method steps of: 

smoothing (1 705) the composite baseline to form an intermediate baseline; and 
forming ( 1 708) a filtered signal (1801) having a second noise component attenuated 

relative to the first noise cnmnonfnt hv <nihrrartino rk^ rnmnncitp (u^ImA ->^A *aa\ 

intermediate baseline to, the input signal. 



9. The method for analyzing an input signal as claimed in claim 8, further comprising the 
method step of replacing (1707) the approximate baseline with the intermediate baseline, and 
wherein the method steps of subtracting the approximate baseline, calculating first and second 
threshold values, substituting the input signal into the approximate baseline, smoothing the 
composite baseline, and replacing the approximate baseline are performed iteratively (1706), 
whereby the second noise component is further attenuated with each iteration. 

10. The method for analyzing an input signal as claimed in claim 1, further comprising the 
method step of: 

smoothing (705) the composite baseline (901) to form an improved baseline (1001) 
having a second error component smaller than the first error component 
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1 1 The method for analyzmg an input stgnaJ as claimed in claim 1 0, further composing the 
method step of replaczng the approximate baseline (1002) with the unproved baseline (1001). 
wherein the method steps of subtracting (702) the approximate baseline (1002), calculating (702) 
first (802) and second (803) threshold values, substituting (704) the input s.gnai (,01) into the 
approximate baseline (,002), smoothing (705) the composite baseline (901), and replacing the 
approximate baseline (1002) are performed iteratively (706), and whexeby the second error 
component is further attenuated with each iteration. 



12. The method for analyzing an input sienai as claims in , - 

*> "«*cm uk input signal is 

chromatographic data. 



13 The method for analyzing an input signal as claimed in Cairn 1 0, further comprising the 
method step of integrating (707) the input signal (10I) ^g the improved baseline (1001). 

14 The method for analyzing an input signal as c.atmed in claim 10, .herein the approximate 
baseline (1002) further comprises a data suture b(x,, the input signal (101) (tinner comprises a 
data sttucture s(x). the baseline-adjusted signa, (801 , funher comprises a data stiucrure a<x). and 
wherein the method step of substituting (704) the input signal (1 0 „ for the approximate baseline 
( 1 002, further compnses substituting (704, the value of «x) into b(x, a. aU values of X where 
a(x) falls within the noise band. 



15. A 

program (6?f^ ; r ^^ ]r . v ^^ t v , v 



program storage medium (623) readable by a computer (622), tangibly embodying 
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the accuracy of an approximate baseline (1002) for an input signal (101), the approximate 
baseline (1002) having a first error component, and the input signal having a background level, 
comprising the method steps of: 

subtracting (702) the approximate baseline (1002) from the input signal (101) to form a 
baseline-adjusted signal (801) having base line -adjusted noise; 

calculating (703) first (802) and second (803) threshold values surrounding the range of 
the baseline-adjusted noise and defining a noise band between the first and second threshold 
values; 

substituting (704) the input signal (!0I) for the approximate baseline (1002) wherever the 
baseline-adjusted signal (801) falls within (902) the noise band to form a composite baseline 
(901); and 

smoothing (705) the composite baseline (901) to form an improved baseline (1001) 
having a second error component smaller than the first error component. 

16. A program storage medium (623) readable by a computer (622), tangibly embodying a 
program (626) of instructions executable by the computer to perform method steps for filtering 
an input signal (503) having a background level, comprising the method steps of: 

smoothing (1201) the input signal (503) to form an approximate baseline; 

subtracting (1202) the approximate baseline from the input signal (503) to form a 
baseline-adjusted signal having baseline-adjusted noise; 

calculating (1203) first and second threshold values surrounding the range of the 
baseline-adjusted noise and defining a noise band between the first and second threshold values; 
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substituting (1204) the input signal (503) for the approximate baseline wherever the 
baseline-adjusted signal falls within the noise band to form a composite baseline. 

17. The program storage medium (623) readable by a computer (622) as claimed ,n claim 16, 
wherein the background level includes a first drift component, further comprising the method 
steps of: 

smoothing (1205) the composite baseline to form an improved baseline; and 
subtracting (1206) the improved baseline from the input signal (503) to form a filtered 
signal (1301) having a second drift component attenuated relive to the first drift component 

18. The program storage medium (623) readable by a computer (622) as claimed in claim 16, 
wherein the background level includes a first drift component and a first noise component, 
further comprising the method step of subtracting (1407) the composite baseline from the input 
signal (503) to form a filtered signal (1601 ) having a second drift component attenuated relative 
to the first drift component, and a second noise component attenuated relative to the first noise 
component. 

19. The program storage medium (623) readable by a computer (622) as claimed in claim 16, 
wherein the background level includes a first noise component, further comprising the method 

steps of: 

smoothing (1705) the composite baseline to form an intermediate baseline; and 




33 

forming (1708) a filtered signal (1801) having a second noise component attenuated 
relative to the first noise component by subtracting the composite baseline from, and adding the 
intermediaie baseline to, the input signal (503). 

20. A data analysis system, comprising: 

an input (621) for receiving an input signal (618) having at least one peak and further 
having a background level including background noise; and 

a processor (622) in communication with the input (62 1 ) which 

subtracts (702) an approximate baseline from the input signal to form a baseline 
adjusted signal having baseline-adjusted noise, 

calculates (703) first and second threshold values surrounding the range of the 
baseline-adjusted noise thus defining a noise band between the first and second threshold values, 
and 

substitutes (704) the input signal for the approximate baseline wherever the 
baseline-adjusted signal falls within the noise band to form a composite baseline. 

2 1 . The data analysis system of claim 20, further comprising a chromatograph (6 1 0) connected 
to the input (62 1 ), the chromatograph (6 1 0) adapted for receiving a chemical sample (601 ) 
having one or more components, the chromatograph (610) further having a column (614) for 
separating the chemical sample (601 ) into the components and a detector (617) responsive to the 
components for generating the input signal (618). 
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22. The data analysts system of claim 20. further composing an output dev lce (628), wherein the 
processor (622) further 

smoothes (705) the composite baseline (901) to form an improved baseline (1 001), 
integrates (707) the input signal (618) using the unproved baseline (1001), and 
transmits the integrated signal to the output device (628). 



23. A method of signal analysis substantially as herein described with reference to Figs 
6 to 18 of the accompanying drawings. 



-•». A program storage medium suhsta 

~ j — — n-iumic tu r igb 

6 to 18 of the accompanying drawings. 

25. A data analysis system substantially as herein described with reference to Figs 6 to 
18 of the accompanying drawings. 
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